Radial equation for an accretion flow driven by Poynting flux 
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Abstract 

Using a toy model of a two dimensional accretion disk, we discuss further 
the radial equation of the accretion flow dominated by the Poynting flux. As- 
suming the force-free condition is valid around the accretion disk, the relation 
between the fluid angular velocity of the disk, Qjj, and the angular velocity 

of the magnetic field lines, Qf, is also discussed. 
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I. INTRODUCTION 



The model of Poynting flux has been suggested for describing the astrophysical jets U 
|| and discussed also for gamma ray bursts @ |]j an d the evolution of a black 
hole-accretion disk system [|]. 

Recently the effect of Poynting flux on accretion flows has been formulated relativistically 
in the background of Kerr geometry for the rotating balck hole at the center [IU| . A simplified 



model of a two-dimensional accretion disk located on the equatorial plane is assumed as in 
the non-relativistic formulation [|TJ. To see the effects of the Poynting flux transparently it 
is also assumed to be non-viscous and non-radiative, although the realistic accretion disk is 
much more complicated |T . 



In this brief note, we reformulate the two dimensional accretion model suggested in ref. 
TOfl , where the time- and toroidal part of the energy- momentum conservation has been 



discussed to calculate the accretion rate and the Poynting flux. The radial part of the 
energy-momentum conservation is discussed in detail to see how the angular velocity of 
magnetic field, Qp, is related to that of the accretion flow, Qp. 

The structure of the electromagnetic field with an accretion disk which lies on 9 = n/2 
plane can be reproduced by assigning the surface charge and current on the plane in addition 
to the bulk charge and current distributions outside the disk. 



We can define a conserved current as 12 



J* = J»e(n/2 -9) + r_Q{9 - tt/2) + f, (1) 

where J± is the bulk current densities( subscripts + and — stand for the upper and lower 
half space respectively) which appear in the Maxwell equation, 

F£% = AnJZ, (2) 

which terminates on the disk, and j M is the surface current density defined on the two 
dimensional accretion disk. Using the current conservation, = we can identify the 
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surface current on the disk as [] 

f = ±-[F 9 »]6(6-n/2). (3) 

Following the standard procedure the surface charge density, a e , and surface current 
density, K, on the disk can be obtained by 

K = ±-[B}x§, (5) 

which are the Gauss' law and the Ampere's law on the surface of the accretion disk. The 
discontinuity of the fields at the disk plane due to the surface charge and current leads to 
the Dirac-delta singularity given by 

3^ = 47^5(0-71-/2), (6) 
d e B =AnK x 5(0-7r/2). (7) 

on the accretion disk. 



II. CONSERVATION OF ENERGY - MOMENTUM TENSOR 

The conservation of the energy momentum tensor given by 

= (8) 

has four equations, one for each v . For the steady and axial symmetric case, there are 
two Killing vectors and therefore correspondingly the energy and the angular momen- 
tum(azimuthal component) of the system are conserved. The relevant combinations of 



1 [A] = A + — A_. The electric and magnetic field are defined for ZAMO(Zero Angular Momentum 
Observer) |i~3fl . 
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equations of v = 0, are responsible for the energy and the angular momentum con- 
servation. The radial component of the equations, v = r, determines the orbital motion of 
the fluid. 



As in ref. ||T0[ we simplify the accretion disk to be non-viscous, cool and non-radiative 
with p = T = for the purpose of investigating the effect of the magnetic field transparently. 
Then the stress-energy tensor of the matter is given simply by 

T% = Pm u»u v . (9) 

It is also assumed that there is no mass flow in the direction perpendicular to the disk: 
u e = 0. 

A. Energy and angular momentum conservation 

The Killing vectors, £ M in t — direction and n in — direction 

£" = (1,0,0,0), ^ = (0,0,0,1) (10) 
define the energy flux and angular momentum flux respectively 

£P = -T«»£ V , C^ = -T^ V „, (11) 

which satisfy 

£% = °> £(W = °- ( 12 ) 

It can be decomposed into two parts: contributions from matter (m) and electromagnetic 
(EM) parts: 

° — °m + c EMi *-(0) — *^n(0) + ^EMty) ^ i0 / 

Since we are interested in a steady and axisymmetric case, the derivatives with respect 
to t and vanishes and hence r— and 9— components of the flux are only considered: 



4 



v~9 



-g 

where 

£ m = -pmU u r , E 6 m = (16) 

&em = -^ — {aE § + / W)^* 4m = ~{-aEf + /3wB>)B*, (17) 
47r p 47rp 

1 WA ^ f fl 1 W 

A D ° 5 '-'EMU) ~ ~~ A 

47T p yvj 4tt p 



^em(4>) = --r-—— B ^ B ^ ^EMU) = - ——B^B 6 - (19) 



Since we assume an idealized thin disk, we do not expect any substantial radial-flow 
of electromagnetic energy and angular momentum compared to the flow of the accreting 
matter which has a delta-function type singularity, for example, 

£ r m = - Pm u u r = -^ Uo u r 6{e - tt/2) (20) 

However the radial flow d r (y/—g£g M ) has no delta-type singularity such that a radial flow 
can be ignored in this thin disk approximation. 

On the other hand the flux in the 9 direction has a singularity on the disk due to the 
discontinuity of the electromagnetic fields. Keeping singular parts only , we obtain 

d r (u )M + + 2nrK r (-aE f + (3zuB § ) = 0, (21) 
d r ( U( p)M + + 27rrK r wB § = 0, (22) 

for energy and angular momentum conservation respectively. The mass accretion rate de- 
rived from the rest-mass conservation is given by [|TT] 



M+ = -2nra m u r (23) 

Eq. (pip and ( |22|) have been derived from the condition that there should be no accumulation 
of energy and angular momentum in the stationary accretion flow driven by Poynting flux 
as discussed in ref. . 



B. Radial equation 



The radial component of the energy-momentum conservation 

= (24) 
is useful to investigate the angular velocity of the fluid. The matter part can be written as 

T£ .,» = - P rn9 rr u%(u,)(^ + ^) (25) 
And the electromagnetic part of the radial equation is given by 

Km ; , = ^d^v^gr&i) + Vem- ( 26 ) 

Since only singular parts are relevant for the thin accretion disk, we get 

1 ~ Sr 



T^M .„ - (27) 

A 'a e E f -K*B B )6{e-Tc/2), (28) 



where — > indicates the procedure of keeping singular parts only. 

After a short straight forward calculations, the radial equation, eq.(23), can be written 

by 



- ^g"u d r (u^(Q D + f^) - ^(v e E f - K*B 8 ) = (29) 

where Qd = u^/u° is the angular velocity of the fluid in an accretion disk. The last term 
corresponds to the Lorentz force exerted by the surface charge and current. 

In the absence of the external magnetic field, the radial equation can be satisfied with 

O = -p-, (30) 

which is one of the characteristics of Keplerian orbit around the black hole. We can see that 
the orbits of the particle in the disk under external electromagnetic field deviate from the 
Keplerian orbit as expected. Hence, the trial ansatz, Qd — used in the literature as 
a first approximation for the order-of-magnitude calculations, might not be valid for a disk 
dominated mainly by the Poynting flux. 



III. FORCE-FREE MAGNETOSPHERE 



In this section we assume that there is a sufficient ambient plasma around the disk to 



maintain the force- free magnetosphere fl| [14]. The force- free condition for the magneto- 



sphere with the current density J M is given by 

= 0. (31) 
With force-free condition, magnetic field lines are spiraling rigidly to form magnetic 



surfaces |13| with the angular velocity fl F and we get 

E = -v F xB p , v F =-{tt F + /3)j> (32) 



a 



where v F is the toroidal velocity observed by ZAMO. Then the toroidal surface current 
density is given byQ 

K+ = a e /v F (33) 

which seems to look strange when compared to the naive expectation 

~ a e v F . (34) 

It is simply because the surface charge defined by the boundary condition is not the charge 
relevant for the surface current and v F is not the velocity of the fluid in the disk. 
Now the Poynting flux is given by 

4m = ^-f£(</>)5 (35) 

and eq.(pl|) becomes 

d r (u )M + - 2TrrK r wQ F B § = 0. (36) 



2 It is consistent with the non-relativistic treatment, Eq.(3.34) multiplied by a and Eq.(3.36) in 
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When compared with eg. fl22|) , the angular velocity of the magnetic surface in a force-free 
magnetosphere is found to be related to the fluid velocity in the following way 

a, = -p£- (37) 
It is inetersting to note that it has no explicit dependence on the electromagnetic field. 



One of the interesting cases is u and are those of the Keplerian orbit [|Kj where 

duo/dr 



du^/dr 

and therefore 



Q K (38) 



n F = Vt K (39) 



However it is subject to the radial equation, Eq . (|29|) , in which the naive suggestion, Eq. (|39j) , 
is not guaranteed in a simple way as discussed in the previous section. In fact the radial 
equation provides more information rather on the relation between and Qp. 
Eq.(|29|) can be written under the force-free condition as 



_ -^ g rr u o 9rU ^Q D _ q f ) _ ^(a.vp - K*)B e = 0. (40) 

The second part represent the Lorentz force in the radial direction and it should be noted 
that the force-free condition does not apply on particle motions in the disk. Hence one can 
see that the angular velocity of the field lines is not the same as that of the disk. It is quite 
contrary to the conventional assumption, Vtp = Qp, for the highly conducting accretion 
disk. It is rather similar to the case with black hole in a force-free magnetosphere, where 
the angular velocity of the black hole, fin, is different from Q F |TB|. This equation shows 



how flp m the force-free magnetosphere can be determined by fl^- 

Naively one can consider the case where the toroidal current is entirely due to the 
surface charge density, <r e , rotating with the angular velocity of the fluid, VLd ■ Then the 
toroidal current is given by 



K+ = <Tv+,vi = -(n D +l3) (41) 
a 
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where vfj is the local velocity of the matter observed by ZAMO. Then Eq. ( f4~0|) becomes 

{ _^l g rr u o d + ^ a ^ B ey {QD _q f) = (42) 
r r z a 

which seem to show that the angular velocity of the magnetic field line is the same as the 
fluid angular velocity, 

Vt F = Q D . (43) 

However it does not hold every where on the disk but only for the particular section of the 
accretion disk where 

4(r)=4W = c. (44) 

which might not be useful in the realistic situation. Of course it is due to an assumption 
that the surface current density is entirely due to the rotating surface density defined by the 
boundary condition. 

IV. SUMMARY 

In this work, a model of the two-dimensional accretion flow driven by Poynting flux 



discussed in |TIJ is reformulated to discuss the radial equation of the acrreting matter. Due 
to the Lorentz force, the angular velocity of the fluid motion, Qp, deviates from the Keplerian 
angular velocity Qk- With force-free magnetosphere, we get a relation between Qp and fluid 
motion: ftp — — ■ However it is found that there is a difference between Qp and ftp as 
for the case of rotating black hole in a force-free magnetosphere. In practical applications 
with accretion disk, the discussions can be very much simplified by assuming Qp> ~ Qk 
and/or Qp ~ Qp>. Hence it would be very interesting to estimate the deviations from those 
ideal equalities and it remains as future numerical works to do. 

It is also worth noting that the recent numerical work [14] demonstrates the possibility of 
a Poynting flux dominated accretion flow. Although a non-relativistic formulation is used, 
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it suggests a possibility that at least a part of the accretion flow can be driven dominantly 
by the Poynting flux. 

The author would like to thank Hongsu Kim for useful discussions. This work was 
supported by Korea Research Foundation Grant (KRF-2000-015-DP0081). 
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